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A bstract: 

We study the effect of the density dependence of the scalar and the vector part 
of the nucleonic self-energy in Relativistic Quantum Molecular Dynamics (RQMD) 



on observables like the transversal flow and the rapidity distribution. The stabil- 
ity of nuclei in RQMD is greatly improved if the density dependence is included in 
the self-energies compared to a calculation assuming always saturation density of 
nuclear matter. Different approaches are studied: The main results are calculated 
with self-energies extracted from a Dirac-Briickner-Hartree-Fock G-matrix of a one 
boson exchange model, i.e. the Bonn potential. These results are compared with 
those obtained by a generalization of static Skyrme force, with calculations in the 
simple linear Walecka model and results of the Briickner-Hartree-Fock G-matrix of 
the Reid soft core potential. The transversal flow is very sensitive to these different 
approaches. A comparison with the data is given. 
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1. Introduction 



The dynamics of heavy-ion reactions depends on the mean field and the description 
of binary collisions. Hence, the mean field and the cross sections for binary collisions 
are the main ingredients of transport models like Quantum Molecular Dynamics 
(QMD) m, U and the BUU-type models @, |, §, §, ^ § which are proven to describe 
successfully the dynamics of heavy-ion collisions. Finally, both ingredients should 
be determined from the same parameter free microscopic force. The strong repulsive 
short range correlations of these forces necessitate a non-pertubative self-consistent 
treatment as done in Bruckner theory. 

In the non-relativistic regime such a self-consistent description of heavy-ion re- 



actions was done in the framework of QMD|| [1JJ. In this approach self-consistent 
potentials as well as in medium cross sections are extracted from non-relativistic G- 
matrix elements which are given by Briickner-Hartree-Fock (BHF) calculations. In 
these calculations the Bethe-Goldstone equation including an angle averaged Pauli 
operator was solved for the situation of two pieces of nuclear matter penetrating 
each other. As a bare nucleon-nucleon interaction the Reid-soft-core potential is 
used in this approach. 

In the covariant extension of QMD, i.e. Relativistic QMD (RQMD) [|TJ, [12], 



13| , up to now only static Skyrme forces combined with the so-called Cugnon 
parametrization of the nucleon-nucleon cross section were used. These general- 
ized Skyrme forces are treated as scalar potentials in the framework of Constrained 
Hamilton Dynamics. But the full Lorentz-structure of the nucleon-nucleon interac- 
tion contains large scalar and vector components. Effects of these relativistic forces 
were already studied in the framework of covariant generalizations of models of the 
BUU-type |14| , 15] , |16| , mostly in the framework of the a — u model, but not yet in 



RQMD. 



In order to determine the scalar and vector components of the self-energy of the 
nucleon from realistic forces in a covariant fashion the solution of the Bethe- Salpeter 
equation with Pauli-operator is needed, which can be achieved in a three dimensional 
reduction of the Bethe-Salpeter equation as given by the Thompson equation. Nu- 
clear matter calculations |17|, [L8| in this framework are done by using one boson ex- 
change potentials (OBEP), e.g. the Bonn potential, as the bare nucleon-nucleon in- 
teraction. The self-energy components extracted from these Dirac-Briickner-Hartree- 
Fock (DBHF) calculations show a strong density dependence while the momentum 
dependence is moderate. Therefore, in the present approach we like to study as a 
first step the density dependence of the self-energy components in a local density 
approximation (LDA). For this purpose we use in RQMD simulations self-energy 
components as extracted from G-matrix elements determined by DBHF calcula- 
tions using the Bonn potential as well as a parametrization fl9| footing on BHF 
calculations using the Reid-soft-core potential. In both approaches we determine 
observables which reflect on the time evolution of the phase space, namely, the 
rapidity distribution and the transversal flow inside the reaction plane. 

In the second chapter we sketch our RQMD approach. In addition, in this chapter 
we give shortly an outline how the results extracted from the Bruckner approach for 
the mean field using the solutions of the non-relativistic Bethe-Goldstone equation 
and the relativistic Bethe-Salpeter equation as well are used within RQMD in a local 
density approximation. In the third chapter we discuss the results and compare them 
with different approaches and data. In the last chapter we finally summarize our 
main conclusion. 
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2. Density Dependent Interactions within RQMD 

In this section we give a brief description of the formalism of relativistic Quantum 
Molecular Dynamics (RQMD) with relativistic self-energies. 

In the formalism of Constrained Hamiltonian Dynamics |2U], [H]], the original 8N 
dimensional phase space of N interacting relativistic particles is reduced to usual 
6N dimensions by 2N constraints. These constraints fix the individual energies and 
parametrize the world lines by fixing the relative time coordinates. 

The energies are fixed by N on-shell conditions 

Ki = (pf - - - (rm - £h 2 = 0, (1) 

which contain scalar and vector potentials. Hence, the constraints (1) take the 
full Lorentz structure of the NN interaction into account explicitly. We regard 
the self-energy components E,- and Ef as sums of two-body mutual interactions- 
at-a-distance, i.e. E^ = J2j^i^\f and Ef = X^jE^. In order to use the self- 
energies provided by relativistic NN interactions, we adopt a kind of local density 
approximation (LDA). In the relativistic Hartree approximation the self-energies are 
proportional to the scalar density p s and the baryon current B^, i.e. 

£(*) = T s p s , W = T V B^. (2) 

Interpreting the sum of the Lorentz scalar two-body interaction densities 

PH = (4^)3/2 ex P(gV4lO, (3) 

as a scalar density p\ = J^j^i Pij one can define in addition the baryon current in a 
consistent way and gets the self-energies of a particle i as 

s, w = r.E^. E? = r B 5>,rtf. (4) 

Here = (g^ - ptjPij)^ Qij = <£ ~ Pij = Pi + Pj , Pij = Pij I (pfj) 1/2 and 
being the 4- velocity of particle j. In this sense the self-energies are given by sums of 
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two-body mutual interactions-at-a-distance. The dependence of these self-energies 
on 5^ ■ ensures the multi-time description used in RQMD without further additional 
approximations. A detailed discussion proving the reliability of these definitions is 
given in the appendix. 

In the case of realistic NN interactions the strong short-range components re- 
quire to take account for the effects of two-body correlations. This can be attained 
in the framework of Bruckner theory. Results of self-consistent Dirac-Briickner cal- 
culations can be used in the approach discussed above with help of effective cou- 
pling functions T s and r„, which are density dependent. This density dependence 
is usually not included. We shall show that the inclusion of this density depen- 
dence improves the description greatly. In this way, realistic forces extracted from 
Dirac-Briickner calculations can be implemented in RQMD. 

For the actual simulations we have used 4 kinds of self-energies extracted from 
the following models and theories: 



(1) Skyrme model for the soft equation of state |L2], [L3| (Skyrme), 

(2) a-uj model H, |3f (QHDI), 



(3) Bruckner Hartree Fock theory with Reid soft core potential |19] (BHF), 

(4) Dirac-Briickner Hartree Fock theory with Bonn potential p4| (DBHF). 

In order to obtain fully covariant constraints, we have fitted the density de- 
pendent coupling constants as functions of the local scalar density to reproduce the 
original calculations of scalar and vector self-energies of models (3) and (4) by means 
of the Lagrange interpolation method. The quality of the fit is shown in Fig.l. 



We have followed the ref. 131 for the choice of the other additional N constraints 



to specify the world lines of the particles: 



Xi = w ijPij1m = °> * = !> -> N ~ 1 ( 5 ) 
Xn = P^Q, - r = (6) 
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with = P^/{P 2 f/ 2 , P^ = EiPi, = EiQi/N and the dimensionless scalar 
weight function 

w ij = -2TT ex P( ( li j / L c) (7) 

Qij/^c 

with L c = 8.66 fm 2 and r denotes a global time evolution parameter. 
The total Hamiltonian is now expressed as 

2JV-1 

H= E ^ (8) 

i=i 

where fa = K { for % < N and 0; = Xi-N for N + 1 < i < 2N — 1, then the equations 
of motions are given by 

dq^/dr = [H, q? ] , dpf/dr = [H, pf] (9) 

where [A, B] means the Poisson bracket of phase space functions A and B. The evo- 
lution of the system can be computed by integrating the set of above equations (9), 
but we have to determine the unknown Lagrange multipliers Aj(r) to calculate the 
dynamics. This can be done using the fact that the complete set of 2N constraints 
(pi ((p2N = Xn) must be fulfilled during the whole time evolution, i.e. 

dfa/dr = dfa/dr + [H, fail = 0. (10) 

Using the Hamiltonian (8), one gets for the Poisson bracket 

2JV-1 

[H,fa}= £ Xj[fa, fa] (11) 
i=i 

thus inverting the matrix Cji = [fa, fa], we obtain the multipliers as 

A, = C- N \ (12) 

If the Dirac's first class condition is fulfilled, i.e. [K i: Kj\ = 0, (in our case the 
condition has been confirmed numerically during the whole time evolution.) the 
Hamiltonian is reduced into a simpler form: 

N 

H = J2 X i R i ( 13 ) 

1=1 



A, = (14) 

where A^ = [Kj, Xi] is a submatrix of the constraint matrix C^. 

For more details, especially for the foundations of these formalisms, we refer to 
HI in case of QMD and to JTIJ and in case of RQMD. We note that the binary 



collisions are dealt in RQMD in the same way as in the non-covariant QMD by using 
Monte Carlo methods. 

3. Results 

First we compare the stability of RQMD nuclei described with different interac- 
tions obtained from the above mentioned 4 models and theories. It is important to 
check the stability at actual boosting energies, since the momentum dependence of 
the self-enegies hides the trivial analytic expressions of covariance. We adopt the en- 
ergy at 1.8 GeV/A for later simulations. Figure 2 shows the root mean square radii of 
Ca versus time, with 6 different approaches: Figure 2 contains the time dependence 
of phenomenological models, i.e. static Skyrme and linear Walecka models. As is 
easily seen, although the Skyrme model works well, but in the simple linear Walecka 
model the nuclei dissassemble after several time steps. This is understandable: the 
strong constant couplings directly reflect on strong repulsive components of the force 
when the local density reaches high values because of the Fermi motion of the nucle- 
ons, however it might be expected that the density dependence of the self-energies 
recovers the discrepancy, since the self-energies are softened in high density domain 
if the two body correlations are included in Bruckner theory (see Fig.l), thus the 
correlation, i.e. density dependence of the couplings would give a counter balance 
automatically inside of nuclei. In order to visualize this mechanism clearly, thus we 
plot the time dependence with microscopic density dependent self-energies and the 
self-energies at the saturation point from BHF and DBHF calculations in Figs. 2. 
One can see that the density dependence of the self-energies strongly improved the 
stability of the nuclei. 



One of the most interesting physical observable in heavy ion collision at inter- 
mediate energies are rapidity distributions and the collective transverse flow. These 
quantities are very sensitive to the interactions used in the models, thus frequently 
used to investigate the nuclear equation of state far from the saturation point. In 
order to be able to compare with the experimental data, we have performed simula- 
tions of the nearly symmetric collision Ar+KCl at 1.8 GeV/A which were measured 
with the streamer chamber at the BEVALAC at LBL p5| . In the calculation we are 
using the above 3 kinds of stable nuclei, i.e. simulated by (1) Skyrme (2) BHF (3) 
DBHF self-energies and compared with the experimental data ||26|| . The initializa- 
tion of the target and the projectile has been carried out by standard Monte Carlo 
methods as described in ref. |T|. According to the ref.^B], we have chosen the impact 
parameter range to be < b < 2.3 fm. The standard Cugnon parametrization is 



used for the binary nucleon-nucleon cross sections |27| . 

The collective transverse flow obtained from our simulation is plotted in Fig.3. 
It is quite remarkable that the conventional Skyrme force (soft EOS) is not able 
to reach the experimentally measured magnitude of the flow distribution at this 
energy. The reason is clear: The force contains only static interactions between 
constituent particles and no explicit energy dependence. However, on the contrary, 
the experimental data lie in between DBHF and BHF results. The Schrodinger 
equivalent potential extracted from the DBHF self-energies shows a linear energy 
dependence in our approach. Hence one gets a repulsive interaction at high relative 
momenta which produce a strong sidewards flow, which is reflected on the P x /A. 
There is still deviation between the theory and the experiment, it might be overcome 
if one takes the configuration dependence (deviation of the momentum distribution 
from a Fermi sphere) in the momentum space of the selfenergies into account, means 
by using local configurations which can be parametrized by two penetrating fermi 
spheres ||19 . 



8 



The rapidity distributions is also shown in Fig.4. The asymmetric property of 
the distribution is due to the system we study. One can notice that the Skyrme 
result exhibits less and BHF gives the highest stopping, but the differences between 
the various approaches are rather small. We found that the Skyrme provides a 
little bit less Pauli blocking (around 10 percent compared with the realistic forces.). 
And also, the stability of single nuclei might affect on this quantity, since it is not 
quite stable using BHF in comparison when using static Skyrme interactions or 
interactions extracted from DBHF (see Fig. 2). 

4. Conclusions 

We investigated the effect of realistic forces in comparison to phenomenological 
forces on the sidewards flow and the rapidity distribution created during heavy-ion 
reactions by using the covariant RQMD approach. The realistic forces are extracted 
from Dirac-Briickner calculations and non-relativistic Bruckner calculations as well 
while the phenomenological forces are determined by local static Skyrme interactions 
or in the framework of a simple linear Walecka model. 

It was proven that single nuclei are stable during the time span needed for a 
simulation of heavy-ion collisions if the full density dependence (beyond a linear 
one) of the interaction is included when using realistic forces. 

The rapidity distribution at the final state shows that realistic interactions create 
slightly more stopping than the static Skyrme interaction. 

Due to the energy dependence and the correct density dependence of these real- 
istic forces the flow is enhanced in comparison to calculations using static Skyrme 
interactions. Furthermore, the results obtained using realistic forces are much closer 
to the data. The agreement with the data might be improved if the self-energies are 
determined including the full local momentum distributions (going beyond a Fermi 
sphere) as it will be done in the near future. 
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Appendix The usual density used in a non-relativistic framework is denned 
as a sum over N Gaussians describing the N nucleons [Q, i.e. 

1 N 

L = 1.08 fm 2 fixes the width of a Gaussian. Reasonable values for p G within the 
multi-time description of RQMD are achieved by propagating all nucleons (forward 
or backward) to the same time coordinate in the frame in which p G has to be 
determined, e.g. in the CMS. 

In order to prove the reliability of the definitions for the scalar density ps and 
the baryon density pB used in eqs.(4) we determined the density profile of pc,Ps 
and ps of a Ca nucleus in its rest frame and for a boosted situation (boosted to 
the CMS energy of a symmetric collision at 2 GeV/A) as well. The results of these 
test calculations are shown in Fig. 5. In the rest frame of the nucleus pc and ps 
coincide as it should be while ps is smaller in the central region. In the boosted 
situation p G and ps are enhanced by a 7— factor compared to ps- Pg an d Pb show a 
quite similar shape in the boosted situation which justifies that the definition used 
to determine ps in RQMD simulations is reliable, pq shows not the full Lorentz con- 
traction because its determination in RQMD simulations contains approximations 
as mentioned above. The reader should note that the density pa given at a special 
point f has half the width of the so-called interaction density evaluated at the center 
of a Gaussian r*j of a particle i. Hence, ps and ps as well are also determined by 
using half the width of the corresponding interaction density when evaluated at a 
special point in space time as done in the calculations shown in Fig.5. 
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Figure Captions 



Fig. 1: Effective Couplings as functions of the scalar density. Solid lines asign 
DBHF(Dirac Brueckner Hartree Fock)/BHF(Brueckner Hartree Fock) calculations, 
dashed lines show the fit used in RQMD simulations. QHDI denotes the values from 
the a — to model at the saturation point. 

Fig. 2: Root mean square radii of single Ca nuclei boosted to the CMS energy 
of a symmetric collision at 1.8 GeV/A as a function of time. The different ap- 
proaches DBHF, BHF and QHDI are explained in the caption of figure 1. 'sat.' 
indicates for T s (p s ) and T v (p s ) at the saturation density. 

Fig. 3: Flow distribution versus rapidity at the final state (60fm/c) of Ar+KCl 
collision at 1.8 GeV/A and b < 2.3 fm. The data are taken from pfi"[ . 

Fig. 4: Rapidity distribution at the final state (60 fm/c) of Ar+KCl collisions 
at 1.8 GeV/A and b < 2.3 fm. 

Fig. 5: Usual density pa, scalar density ps and baryon density ps (definitions 
see text) for a single Ca nucleus in the rest frame of the nucleus (left part) and in 
the CMS (right part). 
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